Ionization energy based Fermi-Dirac statistics 
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Quantitative differences of Lagrange multipliers between standard Fermi-Dirac statistics (FDS) 
and Ionization energy {Ei) based FDS (iFDS) are analyzed in detail. It is shown here that iFDS 
is degenerate and its total energy remains the same with the standard FDS. The total energy can 
be obtained by recasting the Ei in quantized form, as required by the new restrictive condition. 
Basically, iFDS provides an alternative route to calculate the Fermions' distribution spectrum. 
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1. Introduction 

iFDS has been used recently to understand the evo- 
lution of resistivity curves with doping and temperature 
fr) in a wide variety of electronic matter M S S 
0, la 13 ■ Such understanding and iPDS's prediction 
can only be captured by introducing the parameter, Ej. 
Here, the Lagrange multipliers, degeneracy and the total 
energy requirement associated with Ej in iFDS is dis- 
cussed in detail. Both FDS and IFDS are for the half- 
integral spin particles such as electrons and holes. Its to- 
tal wave function, ^' has to be antisymmetric in order to 
satisfy quantum- mechanical symmetry requirement. Un- 
der such condition, interchange of any 2 particles {A and 
B) of different states, Vi and [j ^ i) will result in 
change of sign, hence the wave function for Fermions is 
in the form of 



The negative sign in Eq. that fulfils antisymmetric 
requirement is actually due to one of the eigenvalue of 
exchange operator P = —1. The other eigenvalue, P 
= -1-1 is for Bosons. Ca and Cb denote all the neces- 
sary cartesian coordinates of the particles A and B re- 
spectively. Equation ^ is nothing but Pauli's exclusion 
principle. The one-particle energies Ei, E2, E„i for 
the corresponding one-particle quantum states gi, 52, 
Qm can be rewritten as {E^s ± Ei)i, [E^s ± -B/)2, 
{Eis ± Ei)m- Note here that Eis = Einuiai state- It is 



also important to realize that Eis + Ej = E^ 



lectr 



and 



Eis - Ei = Ehoies- Subsequently, the latter {Eis ± Ei)i 
version where i = 1, 2, m with Ei as an additional 
inclusion will be used to derive iFDS and its Lagrange 
multipliers. This ±£'7 is inserted carefully to justify that 
an electron to occupy a higher state N from initial state 
M is more probable than from initial state L if condition 
Ei{M) < Ei{L) at certain T is satisfied. As for a hole 
to occupy a lower state M from initial state N is more 
probable than to occupy state L if the same condition 
above is satisfied. Eis is the energy of a particle in a 



given system at a certain initial state and ranges from 
-I- 00 to for electrons and to —00 for holes. In con- 
trast, standard FDS only requires Ei {i = 1, 2, m) as 
the energy of a particle at a certain state. 



2. Theoretical details 

Denoting n as the total number of particles with ni 
particles with energy (Eis i rt2 particles with en- 

ergy {Eis ± Ej)2 and so on implies that n = ni + n2 
+ ... + Tim- As a consequence, the number of ways for 
qi quantum states to be arranged among ni particles is 
given as 



P{nuqi) 
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nil{qi - ni)\ 



(2) 



Now it is easy to enumerate the total number of ways 
for q quantum states {q = qi + q2 + ■■■ + (Zm) to be 
arranged among n particles, which is 



P{n,q) = l[ 



^ ni\{qi - ni)\ 



(3) 



The most probable configuration at certain T can be 
obtained by maximizing P(n, q) subject to the restrictive 
conditions 



rii — n, drii — 0. 



(4) 



E 



{E,s ± Ei),n, = E,J2iErs ± Ei),dn, = 0. (5) 



The method of Lagrange multipliers |9| can be em- 
ployed to maximize Eq. (j3Jl- Hence, a new function, 
F{xi,X2, ...^J,, X, ...) = / + /^/i + A/2 +... is introduced 
and all its derivatives are set to zero 
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OF dF dF 

^-0^ aX-o- 



(6) 



As such, one can let the new function in the form of 



F = \nP + fi^dn^ + \^{E,s±Ei)dn,. (7) 



After applying Stirling's approximation, dF/drii can 
be written as 



dF 



= \n{qi - Til) - Inn; + /i + X{Eis ± Ei)i 



0. 



(8) 



Thus, the Fermi-Dirac statistics based on ionization 
energy is simply given by 



1 



q.i exp[^j, + \{Eis ± Ei)i] + 1' 



(9) 



Importantly, the total energy, E in iFDS can be ob- 
tained from Eq. lO, which is 



E ^Y.^E,, ± Ei),n, 

i 

i 

= |^[k?,±k?]^|lk^ (10) 
im Am 

k/ = \i.ionized state, and the ± sign is solely to indi- 
cate that the energy corresponds to electrons is — > -\-oo 
while — oo is for the holes, which satisfy the particle- 
hole symmetry. Consequently, Eq. also implies that 
iFDS does not violate the degeneracy requirements. By 
utihzing Eq. © and taking exp[/^ -|- \{E ± Ei)] ^ 1, one 
can arrive at the probability function for electrons in an 
explicit form as 



feCkis) = exp 



— /i — A 



21.2 



2m 



E, 



(11) 



Similarly, the probability function for the holes is given 



by 



fhi^is) = exp 



Ai + A 



2m 



E, 



(12) 



The parameters and A are the Lagrange multipliers. 
h = h/2'K, h = Planck constant and m is the charge 



carriers' mass. Note that E has been substituted with 
h^'k'^/2m. In the standard FDS, Eqs. ^ and (O are 
simply given by, /e(k) = eyip[-^-X{h^\s? /2m)] and //^(k) 
— exp[/i 4- A(fi,^k^/2m)]. Equation (0J can be rewritten 
by employing the 3D density of states' (DOS) derivative, 
dn = \/k|^dkis/27r^, Eqs. JTTl) and (O, that eventually 
give 



V 
27r2 

V 



e I k exp 



oc 



- A 



2m 



dk 



2n 



2^ ^ I k^exp 



V 
2^2 



^e-^~'^' / ki exp 



A^-A^ 
2m 2m 

h'kl 



dk, 



A- 



2m 



V 

2^' 







j k^exp 




2m 



rfk,,, (13) 



dk,s. (14) 



The respective solutions for Eqs. ifT^ and ((Tl|l are 

/2^A;i2x3/2 



fi + XEi 



ji-XEi 



In 



In 



V \ m 



V\m 



(15) 



(16) 



Note that Eqs. (|15() and H16|) simply imply that 
lie{iFDS) = ii{T = 0) + XEi and Hh{iFDS) = fi{T = 0) 
— XEj. In fact, fi{FDS) need to be varied accordingly 
with doping, on the other hand, iFDS captures the same 
variation due to doping with XEj in which, fi(T — 0) is 
fixed to be a constant (independent of T and doping). 
Furthermore, using Eq. one can obtain 



E = 





oo 

e /ktexp 



4m7r2 
4m7r2 
4m7r2 



2X 



2i,2n 



ft2k' 



2m 



dk 



2m 







1 k^exp 


2™ 





m 



2m 



dk,.. 



dki. 



2A' 



.-l-i{FDS) 



2TTXh? 

3/2 



3/2 



27rA?i2 



(17) 
(18) 
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Again, Eq. H17|l being equal to Eq. (|18|l enable one 
to surmise that the total energy considered in FDS and 
iPDS is exactly the same. Quantitative comparison be- 
tween Eq. H17|l and with the energy of a 3D ideal gas, E 
= 3nkBT/2, after substituting Eq. ^ into Eq. lO will 
enable one to determine A. It is found that A remains the 
same as l/ksT. 

3. Discussion 

Recall that the Ej here corresponds to the energy 
needed to ionize an atom or ion in such a way that 
the electrons are excited to an energy level distanced 
at r, not oo. However, the proportionality, Ei(r = r) 
(X Ei(r = oo) is valid, which has been used to describe 
the experimental data of strongly correlated matter with 
minuscule substitutional doping P, Q, S IS S 0|- Ba- 
sically, at constant temperature (T > 0), FDS predicts 
the distribution spectrum if E is varied, relying on ex- 
ternal inputs such as band gap (Eg) and/or Fermi level 
{Ep(T)). On the other hand, iFDS needs only Ej as an 
external input to predict the variation of E, without re- 
lying on Eg and/or Ef{T) at all, and subsequently its 
distribution spectrum can be obtained as well. Notice 
that Ep comes into iFDS as Ep = constant, indepen- 
dent of T and doping. Ep denotes the Fermi level at 
K. The Ej is microscopically defined as Q 



spectrum. Hence, based on the accuracy of these input 
parameters, one can choose either FDS or iFDS to be 
used for one's theoretical models. 



4. Conclusions 

In conclusion, the relationship between FDS and iFDS 
in term of Lagrange multipliers has been derived and 
shown clearly. The total energy considered in ioniza- 
tion energy based Fermi-Dirac statistics is as same as 
the FDS. Actually, the total energy has been recast into 
a fundamental form that consists of initial state and ion- 
ized state energies. iFDS's prediction are also remark- 
able in non-free-electron metals, namely High- Tc super- 
conductors Q, H, IE Q IS 1^ ; feromagnets [3 and ferro- 
electrics As such, one has the option whether to 
adopt FDS or iFDS based on reliable input parameters 
in which, Ep{T) and/or Eg corresponds to FDS while 
El is connected with iFDS. 
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e(0,k) = l + -f exp[A*(ii;° -i?/)]. (19) 
k 

e(0,k) is the static dielectric function, while ICg rep- 
resents the Thomas-Fermi screening parameter. A* = 
(127reo/e^)n^rB, eo and n are the permittivity of space 
and principal quantum number respectively while tb de- 
notes the Bohr radius. In fact, iFDS and FDS take dif- 
ferent approach in term of energy levels and Fermions' 
excitations to arrive at the same distribution spectrum. 
In simple words, iFDS is new in a sense that it gives one 
an alternative route to obtain the Fermions' distribution 
spectrum in which, FDS needs Ep{T) and/or Eg while 
iFDS needs only Ei to arrive at the same distribution 
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